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ABSTRACT 

We present the first physical model for the non-spherical intra-cluster gas 
distribution in hydrostatic equilibrium under the gravity of triaxial dark matter 
halos. Adopting the concentric triaxial density profiles of the dark halos with 
constant axis ratios proposed by Jing & Suto (2002), we derive an analytical ex¬ 
pression for the triaxial halo potential on the basis of the perturbation theory, and 
hnd the hydrostatic solutions for the gas density and temperature prohles both in 
isothermal and polytropic equations of state. The resulting iso-potential surfaces 
are well approximated by triaxial ellipsoids with the eccentricities dependent on 
the radial distance. We also hnd a formula for the eccentricity ratio between the 
intra-cluster gas and the underlying dark halo. Our results allow one to deter¬ 
mine the shapes of the underlying dark halos from the observed intra-cluster gas 
through the X-ray and/or the Sunyaev-Zel’dovich effects clusters. 

Subject headings: cosmologydheory - dark matter - galaxies: clusters: general - 
X-ray: galaxies 


1. INTRODUCTION 

While clusters of galaxies are regarded as one of the most important cosmological probes, 
the conventional modeling of the intra-cluster gas is very approximate at best. The most 
popular and traditional description is the spherical isothermal (3 model, which has been 
claimed as a good empirical description for the observed intra-cluster gas, and in fact has 
been widely used in practically almost all statistical analyses of galaxy clusters including 
mass, temperature and luminosity functions, and the Hubble constant (or more precisely 
the angular diameter distance) measurement via the Sunyaev-Zel’dovich (SZ) effect (e.g., 
Birkinshaw 1999, for a review). 


The situation is changing rapidly both observationally and theoretically; recent X-ray 
satellites such as Chandra and XMM-Newton have greatly improved the angular-resolution 
and the statistical reliability of X-ray maps of galaxy clusters, and often detected the signih- 
cant departure from the isothermality and the non-sphericity of the intra-cluster gas. Buote 
et ah (2002), for instance, reported evidence from the Chandra observation for a flattened 
triaxial dark matter halo around the elliptical galaxy NGC720 (see also, Romanowsky & 
Kochanek 1998). In the radio band the bolometer array technique improved the angular 
resolution of the cluster SZ map. One snccessfnl example is the discovery of snbstructure 
of the most Inminous X-ray cluster RX J1347-1145 by Komatsu et al. (2001), which was 
subsequently conhrmed with Chandra (Allen et al. 2002). 

The theoretical nnderstanding of dark matter halos is also advanced signihcantly. The 
spherical averaged density prohles of dark halos are known to be well approximated by a 
seqnence of nniversal density prohles (Navarro, Frenk & White 1996, 1997; Fnkushige & 
Makino 1997; Moore et al. 1998; Jing & Snto 2000). The importance of non-spherical 
effects on the mass fnnction of dark halos was hrst raised by Monaco (1995), and then 
stndied analytically (Andit, Teyssier, & Alimi 1997; Lee & Shandarin 1998), and in detail 
with numerical simulations (Lee & Shandarin 1999; Sheth & Tormen 1999; Jenkins et al. 
2001). More importantly, the systematic effect of the non-sphericity of intra-cluster gas on 
the estimate of the Hnbble constant from the SZ observations is well recognized (Birkinshaw, 
Hnghes, & Arnand 1991; Inagaki, Snginohara & Snto 1995; Yoshikawa, Itoh & Snto 1998; 
Hughes & Birkinshaw 1998; Fox & Pen 2002). In particular the recent work of Jing & 
Snto (2002) showed that the simnlated halos are well approximated by a sequence of the 
concentric triaxial distribntion with their axis directions being fairly aligned. This natnrally 
opens a possibility to compnte the hydrostatic equilibrinm solntion of the intra-cluster gas 
under the gravity of the triaxial dark matter halos. Generalizing the prescription of Makino, 
Sasaki & Snto (1998) and Snto, Sasaki, & Makino (1998) who obtained the hydrostatic 
eqnilibrinm solution for the spherical halo prohles, we hnd a series of analytical solntions for 
gas embedded in the triaxial halo prohles. In the present paper, we focus on the mathematical 
aspects of the problem, and their astrophysical applications (weak/strong lensing, SZ and 
X-ray observations) will be described elsewhere. 

The rest of the paper is organized as follows: In §2 we briehy ontline the procednre to 
compnte gas and temperatnre prohles in hydrostatic eqnilibrinm with gravitational potential 
of dark matter halos. In §3 we present the analytical derivation of the triaxial halo potential 
nsing the pertnrbation theory, and describe the shapes of the iso-potential snrfaces. The 
analytical and nnmerical results for the gas density and temperatnre prohles are snmmarized 
in §4, and §5 is devoted to conclusions and discussion. In Appendix, we provide the detailed 
analytical expressions of our perturbative results. 
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2. DISTRIBUTION OF INTRA-CLUSTER GAS IN HYDROSTATIC 

EQUILIBRIUM 


If the intra-cluster gas embedded in the dark matter halo is in hydrostatic equilibrium, 
its distribution is described by 

ivp, = -V4, (1) 

Pg 

where <h is the gravitational potential of the system, and pg and Pg represent the density and 
the pressure of the intra-cluster gas, respectively. In what follows, we neglect the contribution 
of the gas and stellar masses, and assume that the gravitational potential of the total system 
is well approximated by that of the halo (Suto, Sasaki & Makino 1998). 

Once the equation of state for the intra-cluster gas and the density prohle of the dark 
matter halo are given, then the density and temperature prohles of intra-cluster gas are 
obtained by solving equation ( 1 ). Although the intra-cluster gas is commonly assumed to be 
isothermal, it has been recently claimed that the intra-cluster gas might be better described 
as polytropic (Markevitch et al. 1998; Komatsu & Seljak 2001). Thus we consider both 
isothermal and polytropic cases. 


Consider first the isothermal gas with an equation of state: 

Pg = Kpg, K ^ (2) 

prUp 

where Tg, ks, p, and nip denote the (constant) gas temperature, the Boltzmann constant, 
the mean molecular weight, and the proton mass, respectively. In this case the hydrostatic 
equilibrium equation ( 1 ) is easily solved to yield 


_l($_$o) , 

where <l>o is an integration constant that can be hxed by the condition of Pg{r = 0) 
Note that the value of $0 depends on the value of K. 



( 3 ) 


PgO- 


Turn next to the polytropic model with an equation of state: 

P, = A'opJ, A'o = hlh, (4) 

^ pnip 

where 7 ( 7 ^ 1 ) is the polytropic index, and T^o is the value of the gas temperature at r = 0 . 
In this polytropic model, the gas temperature is not constant but proportional to So, 

in this case, equation ( 1 ) reduces to 


7-1 

^07 


V$, 


V(p'^-') 


( 5 ) 



4 


whose solution is derived as 


TgO 


1-7 

KqJ 




Pl_ 

PgO 


1-7 

Ko'j 


($ - $o) 


l/(7-l) 


( 6 ) 


Thus, the general density and temperature prohles of intra-cluster gas in both isothermal 
and polytropic cases can be straightforwardly obtained in terms of the gravitational potential 
of the underlying halo. Furthermore, equations (3) and (6) also imply that the iso-potential 
surfaces of the triaxial dark halo coincide with the iso-density surfaces of the intra-cluster 
gas. This is simply a direct consequence of X-ray shape theorem (Buote & Canizares 1994; 
Buote et ah 2002); the hydrostatic equilibrium equation (1) yields 


VPg X V<F = Vpg X V<h = 0. 


( 7 ) 


Moreover, with the polytropic equation of state (eq. [4]), the gas temperature also satishes 


VTg X V4> = 0. 


( 8 ) 


Therefore, the gas density, temperature, and pressure are constant on the iso-potential sur¬ 
faces at the same time. In consequence, the gravitational potential is the central quantity for 
the hydrostatic gas model. In the following section we derive an analytical expression for the 
triaxial dark halo potential, and study the general properties of the resulting iso-potential 
surfaces. 


3. GRAVITATIONAL POTENTIAL OF TRIAXIAL DARK HALOS 


3.1. Tricixial Density Profile of Dark Halos 


In what follows, we adopt the following concentric triaxial density prohle for the dark 
matter halos (Jing & Suto 2002): 


p{R) 


hcPcrit 


( 9 ) 


where Rq is the scale length, 6c is the dimensionless characteristic density contrast with 
respect to the critical density pcut of the universe at the present epoch, and a represents 
the inner slope of the density prohle. Note that equation (9) is identical to the one which 
describes the spherical density prohle (Navarro, Frenk & White 1997; Fukushige & Makino 
1997; Moore et al. 1998) except that the spherical radius r is replaced by the the major axis 
length R of the iso-density surfaces : 

Va" 


R^ = 


{a>b>c). 


( 10 ) 
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The N-body simulations by Jing & Suto (2002) imply that the best-£t values of a are 
a = 1 and a = 3/2 for cluster and galactic-scale halos, respectively. Thus we pay particular 
attention to these two cases. 

We quantify the ellipsoidal shape of a halo iso-density surface by defining the two ec¬ 
centricities: 



and a > b > c implies Cb < Oc ■ The values of (a = b, c) measure the degree of the devi¬ 
ation of the ellipsoidal iso-density surfaces from the spherical ones along the corresponding 
principal axis direction. 


3.2. Perturbative Expansion of the Tricixial Halo Potential 


The gravitational potential of a dark halo with the triaxial density profile (eq. [9]) can 
be written as (Binney & Tremaine 1987) 


<h(r) 

'ipim) 


-tiG 


(cx)) — Ipim) 


'o {t + a^){T + b‘^){T + c^) 
p{R)RdR, 


X y z 


dr, 


( 12 ) 

(13) 

(14) 


where r = {x,y,z) = r(sin 6^ cos 0, sin 6* sin 0, cos 6^), and r labels the iso-potential surfaces 
on which m = m(r, r) is constant. While numerical integration is required in general to 
obtain the triaxial gravitational potential with equations (9) to (14), small eccentricities 
(e^ < 1) cases enable us to approach this problem analytically with the perturbative 

expansion. 

The theory for the ellipsoidal perturbation of the spherical systems has a long history 
(e.g., see Jeffreys 1976, and references therein). In performing the ellipsoidal perturbation, 
there are two different ways to arrange the perturbed terms. One is the equal-volume ap¬ 
proach: the perturbed ellipsoids have the same volumes as the unperturbed original spheres. 
The other is the equal-length approach: the perturbed ellipsoids and the unperturbed spheres 
have the same length scales. For instance, one can force the major axis lengths of the per¬ 
turbed ellipsoids to be same as the radii of the original spheres. The two approaches differ 
only in arranging the perturbed ellipsoidal terms, and the perturbation theory is basically 
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the same either way. Though the equal-volume approach is more conventional, we take the 
equal-length approach here, for which the reason is as follows. 

We develop the perturbation theory here mainly to provide an ellipsoidal correction 
to the conventional spherical modeling of the observed clusters. What is observationally 
measurable, however, is not the volumes but the length scales of the clusters. For example, 
in SZ observations, the cluster length scales in the direction perpendicular to the line-of-sight 
are measured. In the standard spherical model, this tangential length scales are taken as 
the spherical radii of the clusters. In ellipsoidal model, however, the same tangential length 
scales are no longer same as the line-of-sight length. Thus, in the two models, the volumes 
of the clusters are not the same. Thus, given the future application of the perturbation 
theory to the observed clusters, we perform the perturbative expansion of the gravitational 
potential in such a way that the perturbed ellipsoids have the same length scales (here, the 
major axis lengths) as the radii of the spheres but do not have the same volumes. Derivation 
of the perturbative expansion of <F(r) up to the hrst order of in this way proceeds as 
follows. 


First, dehned in equation (13) for the given triaxial density prohle p{R) in equation 
(9) is analytically given as 


,, , 2FllS,p„;, , floV" 

Expressing b and c in terms of a as 6 = Oa/I — and c = Oa/I — by equation (11), we 
expand m in equation (14) with respect to and to their its firstorder: 


m ~ TUsp 




el sin^ 9 sin^ 0 + cos^ 9\ 
2 


(16) 


where = ar / ya? + T is the value of m for the spherical case (e^ = Cc = 0). 
Equations (15) and (16) yield the first order approximation of ip{m) : 


^jJ{m) ^ 2Rl5cPcnt 1 + 


Rn 


m. 


a-2 


1 ^ Cfe sin^ 6 sin^ 0 + cos^ 9 


2 — a ' rrisv + Rq 2 


^sp / ^ ••^sp 

Similarly, we approximate \/{t + a‘^){T + b‘^){T + c^) and be/a in equation (12) as 

^ I \ 

_ 

a 


V(T + a^)(T + b‘‘){T + e) ~^(i + ^"" 


^3^3 


a I 1 — 


el + el 


( 17 ) 

( 18 ) 


Now that all the quantities in the integral of equation (12) are expressed in terms of 
nisp as equations (16) to (18), we change the integration variable from r to nisp, and perform 
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the integration of equation ( 12 ) to obtain the following Ist-order approximation for <h: 


$(u) ^ C 


i"i(w) + 


sm^^sin^0 + eg 


u] 


(19) 


where u = r/i?o; and C = 47rG(5cPcrit-Ro) and the three functions, -Fi(m), F 2 (M),and F 2 ,{u) are 
dehned as 


Fi{u) ^ 

F2{u) = 
F^iu) = 


a — 2 
1 

a — 2 


i_i r/ i 


U 


'0 


t + 1 


2-a 


dt 


3 u 


'0 


t + 1 


2-a 


dt -r 


1 r 


(t + 1) 


2-a 


dt 


1 r 




{t + 1 ) 


3—0; 


dt. 


( 20 ) 

( 21 ) 

( 22 ) 


Equation (19) is valid for any arbitrary value of a. For the interesting values of a = 1 and 
a = 3/2, Fi, F 2 and F^ can be written in terms of elementary functions (see Appendix A). 

Here, Fi{u) represents the spherical contribution to *h(u), i.e., CFi{u) = 47rG6,p,,itRlF,{u) 
is nothing but the gravitational potential for the case of the spherical density profile (Makino, 
Sasaki & Suto 1998; Suto, Sasaki & Makino 1998). F 2 represents another spherical contri¬ 
bution that has arisen due to the volume changes of the perturbed ellipsoidal density prohles 
from the spherical ones, while F^ represents the non-spherical deviation of *h(u) from the 
spherical potential. If we take the equal-volume approach in performing the perturbation, 
then we will end up with exactlly same Fi{u) and F^i^u) but no F 2 . Figure 1 plots these 
functions in a wide range of the rescaled radius u for a = 1 and a = 3/2. As shown, the 
magnitude of F^ is an order of magnitude smaller than those of Fi and F 2 . Since the de¬ 
pendence of *h(u) on 6 and </> comes from the Fs-term (eq.[19]), it implies that *h(u) is fairly 
insensitive to 6 and 0. The asymptotic limit of F 3 /F 1 at m 00 has been calculated to 
converge to zero, which implies that the triaxial potential is almost indistinguishable from 
the spherical one at large distance (m 3> 1 ). 

Equation (19) is valid only for e 2 < 1 . Jing & Suto (2002), however, showed that the 
axes ratios of cluster scale halos are typically 0.4 < c/a < b/a < 0.8, which corresponds to 
0.3 < Cb < Cc < 0.8. To examine how well equation (19) works in this realistic range of e„, 
we integrate equation (12) numerically to obtain the real gravitational potential. Figure 2 
compares the perturbative potentials (dashed lines) to the numerical results (solid lines) for 
various values of the halo eccentricities both in a = 1 and a = 3/2 cases. It is clear that 
the perturbative potentials agree with the numerical results excellently. We also calculate 
the ratios of the perturbative potentials to the numerical ones, and hud the ratios less than 
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1.1 even for eb = 0.6 and Cc = 0.8, which indicates that equation (19) is indeed a good 
approximation. 

Incidently, the gravitational potentials for the axis symmetric halos with the single 
eccentricity e can be also obtained from equation (19): 

^(u)^^ F^{u)±j{F2iu) + cos‘^eF^{u)} , (23) 

where the positive and the negative signs in front of e^/2 correspond to the oblate [a = b > c, 
e = a/ 1 — (c/a)^)] and the prolate [a < b = c, e = a/ 1 — (a/&)^] cases, respectively. Note 
that ; 2 -direction is always chosen as the symmetric axis in the above expression. 


3.3. Iso-Potential Surfaces and Potential Profiles 


Strictly speaking, the iso-potential surfaces of ellipsoidal dark halos are not necessarily 
exact ellipsoids (Binney & Tremaine 1987). However, we have found that the iso-potential 
surfaces are still best well approximated as ellipsoids (see also Figs.4 and 5 below). Thus, 
we model the iso-potential surfaces as triaxial ellipsoids with the rescaled major axis length, 
and the two eccentricities, et and ec- 

+ r^) = (24) 

If equation (24) is a good and consistent approximation for the real iso-potential surfaces in 
the frame of our perturbative approach, we should be able to find the potential profile, <F, 
that satishes *h(u) = ‘F(.C). To hnd a functional form of <h, we hrst expand equation (24) 
assuming 1: 

Substituting equation (25) into ‘h(O) and expanding to hrstorder of e^, we have 


$ M 1 -h 


el sin^ 9 sin^ 0 + cos^ 9 


= I>(m) 


(9<F(m) sin^ 6* sin^ 0-f cos^ 6* 
du “ 2 ■ 


Comparing equation (19) with equation (26) indicates that the functional form of <h is 


” 2 2 " 
i(i) = C fiK) + , 


( 27 ) 
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and that the eccentricities of iso-potential snrfaces are written in terms of their halo 
eccentricities as 

4 ^ F^u) 

el uduFiiu) 1 — {a — 2)Fi{u) — ’ 

Note that e^- for the gravitational potential depends on u unlike the constant e^j for the 
adopted dark matter halo prohle. Figure 3 plots eo-/eo- as a function of n for a = 1 and 
a = 3/2 cases. The a = 1 case has a higher ratio of to than the a = 3/2 case. In the 
whole range of u, €„/€„ is less than unity, and decreases mildly as u increases. 


{a - 2)^3 


u 




4. GAS DENSITY AND TEMPERATURE PROFILES 

In §2, we showed that the gravitational potential of dark halos is the central quantity for 
the hydrostatic intra-cluster gas distribution, and that the gas iso-density and temperature 
surfaces coincide with the halo iso-potential surfaces. Thus, one can restate every feature of 
the iso-potential surfaces found in §3 in terms of the iso-density surfaces of the intra-cluster 
gas: the iso-density and temperature surfaces of the intra-cluster gas in triaxial dark halos 
are approximately triaxial ellipsoids whose eccentricities are related to that of the underlying 
halos via equation (28), decreasing with radial distance. 

Note that we have derived equation (28) using three different approximations: the 
hrst order approximation of the gravitational potential with the perturbation expansion 
(eq.[19]); the approximation of the iso-potential surfaces with triaxial ellipsoids (eq. [24]); the 
hrst order approximation of the iso-potentential surface eccentricities with the perturbation 
expansion (eq.[26]). The accuracy of equation (28) depends on the accumulated errors from 
the three approximations. In §3.2, we have already shown that equation (19) is an excellent 
approximation within the 10% error. 

Figure 4 plots the 2D contours of the gas iso-density surfaces (solid lines) and the 
ellipsoidal approximations (dot-dashed lines) on the equatorial plane (6* = 7i/2) at three 
different radii r = 0.1Ro,-Ro and 10i?o for the given halo eccentricities Cb = 0.6 and Cc = 
0.8. The contours of the gas iso-density surfaces are found numerically directly from the 
gravitational potential given as equation (12). We also plot the 2D shapes of the underlying 
halos (dashed lines) on the equatorial plane, and the circles (dotted lines) as well for reference. 
Figure 4 shows that the approximation of the gas density surfaces (or equivalently,the iso¬ 
potential surfaces) as ellipsoids work quite well even for large halo eccentricities e?, = 0.6 and 
Cc = 0.8. It also demonstrates explicitly that the intra-cluster gas is more spherical in the 
outer part than in the inner part of the potential, and that it is overall more spherical than 
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the underlying dark halo. This is intuitively understood because the potential represents 
the overall average of the local density prohle, and also because the gas pressure is isotropic 
unlike the possible anisotropic velocity ellipsoids for collisionless dark matter. This fact has 
been also detected by hydrodynamic simulations (e.g., see Fig. 4 in Yoshikawa, Taruya, & 
Suto 2001). 

Figure 5 plots the same as Figure 4 but for the case of halo eccentricities, e?, = 0.8 and 
Cc = 0.8. It shows clearly that the ellipsoidal approximations (dot-dashed lines) still work 
fairly well but that there are some noticeable difference between the perturbative (solid lines) 
and the numerical (dot-dashed lines) results. To quantify the difference of the ellipsoidal ap¬ 
proximations of eqnation (28) to the gas iso-density snrfaces, we measure the eccentricities 
of the gas iso-density surfaces directly from the contours of the numerically calculated grav¬ 
itational potential, and determine the ratios between these nnmerical eccentricities and 
analytical approximations Fignres 6 and 7 show the results as functions of the halo ec¬ 
centricities. Note that depends on both e?,, and Cc, so that ^ even when c;, = 0 
if Cc 7 ^ 0. In contrast, for the analytical approximations of equation (28), e?, depends only 
on Cfo. Figures 6 and 7 also snggest that is not so sensitive to the values of u and 

a. On the whole, the values of are less than 1.2, indicating that the error involved 

in approximating ej)"™ to is less 20% at most. Furthermore, to the approximation error, 
we provide the following htting formula to 

1 + [0.1 + 0.05 log(l + u)]el + [0.2 + 0.03 log(l + u)]el (29) 

1 + [0.1 + 0.09 log(l + u)]el + [0.2 + 0.03 log(l + u)]el. (30) 

Equations (29) and (30) represent the accuracy of equation (28), qnantifying the calcula¬ 
tion errors accnmulated from all the above three approximations made in the derivation of 
eqnation (28). In summary, equation (28) is accurate within 10% errors for < 0.6, while 
within 20% errors for < 0.8. 

By equations (3),(6) and (19), we obtain the density and temperatnre prohles for both 
the isothermal and the polytropic gases in the triaxial dark halos, and plot the hnal results 
in Figure 8 for a = 1 and a = 3/2 cases. For this hgure, we choose c;, = 0.6, Cc = 0.8, 
7 = 1.15, and (1 — 7 )/(iFo 7 ) = 1. For comparison, we also plot the axis-symmetric (both 
oblate and prolate with e = 0.8) and the spherical cases. 

Both in isothermal and polytropic cases, the significant deviation of the resulting prohles 
from the sphericity are manifest especially in the range of r > Rq, increasing with r. The 
a = 3/2 case shows non-negligible deviations even for r < 10“^i?o- If is interesting to note 
that the oblate and prolate prohle curves are symmetric about the spherical ones This is due 


^num 

_ 

^per 

^num 

per 
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to the fact that the gravitational potentials for the oblate and the prolate halos differ only 
by the sign before the non-spherical perturbative term when the ; 2 -direction is chosen as the 
symmetric axis (see §3.1). 


5. CONCLUSIONS AND DISCUSSION 

We have adopted the triaxial halo density prohle suggested by Jing & Suto (2002), 
and calculated the halo gravitational potential by using perturbation, assuming the small 
eccentricities of the underlying halo. The approximations have been shown to be valid even 
for the halos with fairly large eccentricities. With the resulting halo potentials, we found the 
solutions to the hydrostatic equilibrium gas equations both for the isothermal and polytropic 
gases. The corresponding gas density and temperature prohles have been shown to deviate 
from the conventional spherical models signihcantly at large radial distances. 

We have also derived a useful analytical formula for the eccentricity ratio between the 
gas and the halo as a function of radial distance. It has been shown that the intra-cluster 
gas is rounder than the underlying halo, and that the gas eccentricity decreases with radial 
distance. It provides a quantitative explanation about why the gas inside dark matter halos 
are observed to be less elongated than the halo themselves in hydrodynamic simulations. We 
expect wide applications of this formula: First, it can be applied to clusters to determine the 
shapes of the unseen dark halos from the observed shapes of the intra-cluster gas from the 
X-ray and Sunyaev-Zel’dovich effects. Second, it can be used to reduce the errors involved 
in the measurement of the Hubble constant caused by the cluster asphericity. Third, it can 
be also useful in determining cosmic shear from the weak gravitational lensing surveys, and 
so on. 

It is worth emphasizing that our analytical results are not empirical htting formula, but 
derived from the hrst principles using the perturbation theory. The only assumptions made 
in our derivation is that the intra-cluster gas is in hydrostatic equilibrium. Thus, it will be 
also possible to test the validity of hydrostatic equilibrium of the intra-cluster gas model 
using our analytical results. Furthermore, our results are quite general in the sense that 
they are derived for any arbitrary value of the gas constants, 7 and Kq, and thus can used 
to constrain the values of 7 and Kq through the direct comparison with the observational 
data. 


We thank A. Taruya, M. Oguri, and T. Kuwabara for their useful comments. We 
also thank the anonymous referee for careful reading of the original manuscript and many 
helpful suggestions. J. L. gratefully acknowledges the support from the JSPS fellowship. 
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A. Analytic Expressions for a = 1 and a = 3/2 cases 

For the inner spectral indices of a = 1 and a = 3/2, which are indeed suggested 
from numerical simulations, F) and e„/e„ defined in §3 can be written explicitly in terms of 
elementary functions. We provide their expressions here which may be useful in quantitative 
confrontation of our results with observations. 

For 0 = 1, 


F,{u) 

F2{u) 

Fsiu) 


u 


- ln(l 
u 

1 2u^ 
'3 + 

2 - 3u 


u 


3u + 6 


Qv? 

6 3 


■U" 


+ (-^ + F), 


2u‘^{l + u) 


ln(l u), 


el 6(1 -|- u) ln(l + u) + — 3u^ — 6u 

el 2v?[{l + + u) — u 


(Al) 

(A2) 

(A3) 

(A4) 


For a 

= 3/2, 

Fiiu) 

= 2 + 2^+“ 

V u 

F2{u) 

4/5 5 

“ 3 ^ V4m3 + 12^2 

Fsiu) 

/ 15 5 

\4u3 4^2 : 


el (15 -5u- 


- ln(v/M -h vT+u), 
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(AS) 
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6u 3 
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1 ) \n{y/u + Vl + u), (A6) 
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+ “—n ln(y^vu-j-T), 


1 -I- u Au^ 
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u = r/Ro 

Fig. 1.— The behavior of the three functions, Fi, F 2 and F 3 (solid, dashed, and long-dashed 
lines, respectively) given in the perturbative result of the gravitational potential (eqs.[ 20 ],[ 21 ] 
and [22]). Upper Panel it corresponds to the case where the inner slope of the halo density 
prohle defined in equation (9) has the value of a = 1. Lower Panel it corresponds to 
a = 3/2. 
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u = r/Ro u = r/Ro 


Fig. 2.— Comparison of the numerical gravitational potentials (solid lines) with the pertur¬ 
bative results (dashed lines) on the equatorial plane of 0 = | and 0 = 0. Top\ it correspond 
to the case where the underlying halo has the eccentricities, Cb = 0.1, Cc = 0.3. Middle, 
eh = 0.3, Cc = 0.5. Bottom: eb = 0.6, Cc = 0.8. The left and right panels correspond to a = 1 
and a = 3/2, respectively. 
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Fig. 3.— The ratio of the eccentricity of the iso-potential surface to that of the halo iso¬ 
density surface from the perturbative result (28); a = 1 (solid line) and a = 3/2 (dashed 
line). 
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Fig. 4.— The shapes of the iso-density surfaces of the gas embedded in the halo with the 
eccentricities Ch = 0.6 and Cc = 0.8 at three different radii, r = 0.1i?o,-Ro and 10i?o on the 
equatorial plane {9 = 7r/2). The sizes of the curves are arbitrary. In each panel, the solid 
and the dot-dashed lines represent the the numerical and the perturbative results of the gas 
density prohle, respectively. While the dashed lines represents the halo iso-density surfaces. 
The dotted circles are also plotted for comparison. 
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Fig. 5.— Same as Figure 4, but for larger halo eccentricities,= 0.8 and Cc = 0.8 
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0.5 0.5 


Fig. 6.— The ratio between the gas eccentricity htted to the nnmerical results and the 
perturbative ones as a function of the halo eccentricity e^. In each panel, the solid, the 
dashed and the long-dashed lines represent the cases of Cc = 0.8, 0.5 and 0.3, respectively. 
The dotted line represents the htting formula given in equations (29) and (30). 
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Fig. 7.— Same as Figure 6, but as a function of Cc for the three different values of = 0.1, 0.3 
and 0 . 6 . 
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Fig. 8.— The hydrostatic equilibrium solutions for the intra-cluster gas distribution. Top: 
the density prohles of the isothermal gas. Middle-, the density prohles of the polytropic gas 
with the polytropic index of 7 = 1.15. Bottom-, the temperature prohles of the polytropic 
gas. In each panel, the solid, the dashed, the dotted, and the dot-dashed lines represent the 
gas prohles for the triaxial (e^ = 0.6, Cc = 0.8), the oblate (e = 0.8), the spherical, and the 
prolate (e = 0.8) halo cases, respectively. 















